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We present an heuristic argument for the universal area product: A+A_ = (8nJ) 2 + (47tQ 2 ) 2 for 
a four-dimensional, stationary, axisymmetric, electrically charged black hole with an arbitrary sta¬ 
tionary axisymmetric distribution of external matter (possibly charged), derived by Marcus Ansorg 
and Jorg Hennig. Here A+ and A_ are the areas of the event and Cauchy horizons, and J and 
Q are the angular momentum and electric charge. Based on this argument, we conjecture that a 
universal area product holds for higher-dimensional, stationary, multi-horizon black objects in the 
presence of an external stationary charged distribution of matter. 
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The product of all horizon areas for general rotat¬ 
ing multi-charge black objects in four- and higher¬ 
dimensional asymptotically flat or anti-de Sitter space- 
times depends only on the quantized charges, quantized 
angular momenta, and the cosmological constant [lj. The 
remarkable result of quantization of the product of the 
horizon areas was already implicitly presented in the 
earlier works by Larsen and Cvetic where the en¬ 
tropies of the horizons (proportional to their areas) of 
four- and five-dimensional black holes were expressed in 
terms of the excitation numbers of the left and right mov¬ 
ing modes of a weakly-coupled two-dimensional confor¬ 
mal field theory. The quantization rule of the product 
of the areas follows from this rule Q- Recently Gold¬ 
stein, Jejjala, and Nampuri [§| showed that the prod¬ 
uct of areas for fixed (7(1) charges is also an invariant 
for all static spherically symmetric black holes in four¬ 
dimensional N = 2 supergravity. 

As it was stated by Cvetic, Gibbons, and Pope 0, if 
the cosmological constant is quantized, the area products 
may provide a “looking glass” for probing the microscop¬ 
ies of black objects. The solutions considered in jjj rep¬ 
resent isolated black objects. Thus, it is important to 
study whether the area products survive in the presence 
of external matter and fields. 

Here we present arguments for a conjecture that this 
is likely to be the case. Our conjecture is based on 
a heuristic argument for the universal area product of 
the horizon areas of a four-dimensional stationary and 
electrically charged black hole distorted by arbitrary sta¬ 
tionary axisymmetric electromagnetic and gravitational 
fields which could be due to an electrically charged sta¬ 
tionary distribution of matter around the black hole. 

The relatioifj] A 2 = (8nJ) 2 + (47rQ 2 ) 2 , where A is the 
extremal horizon area, J is the angular momentum, and 
Q is the electric charge, was proven for an electrically 
charged, extremal black hole with arbitrary surrounding 


matter by Ansorg and Pfister [§j. A universal area prod¬ 
uct of the horizon areas of a stationary and axisymmet¬ 
ric black hole surrounded by matter, A + A_ = (8nJ) 2 , 
where A + and A_ are the areas of the event horizon 
(EH) and Cauchy horizon (CH), was constructed by An¬ 
sorg and Hennig in 1 10], and later its validity was numer¬ 
ically illustrated in [llj . Finally, the universal area prod¬ 
uct A + A_ = (87r J) 2 + (47 tQ 2 ) 2 for a four-dimensional, 
stationary, axisymmetric, electrically charged black hole 
with surrounding matter was proved by Ansorg and Hen¬ 
nig 00 . One can check directly that such a relation 
holds for a Kerr-Newman black hole (see, e.g., 0 )- 

It is rather surprising that the same relation holds in 
the presence of distortion of the black hole by, for ex¬ 
ample, surrounding matter. Ansorg and Hennig 00 
proved the validity of this relation by using the inverse 
scattering method applied to the linear matrix problem 
whose integrability conditions are equivalent to two com¬ 
plex Ernst equations representing the Einstein-Maxwell 
system. Here we give an heuristic argument explaining 
such a relation. 

Our model consists of a distorted Kerr-Newman black 
hole of angular momentum J, electric charge Q , and a 
surrounding stationary, axisymmetric, and, in general, 
electrically charged distribution of matter representing 
distortion sources which is defined by a set of parameters, 
dB The space-time beyond the sources is asymptoti¬ 
cally flat, and the black hole is in stationary equilibrium. 

Let us consider a quasi-stationary transition from an 
undistorted Kerr-Newman black hole solution to another 
stationary and axisymmetric but distorted black hole de¬ 
fined by different values of the parameters Di , such that 
the black hole remains close to equilibrium during the 
transition. An undistorted black hole corresponds to van¬ 
ishing Di, e.g. when the distortion sources are located 
at infinity. The quasi-stationary transition consists of 
small transition steps, such that at the end of each step 
the black hole settles down to another stationary state. 
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1 We use a system of units with G = c = 1 and spacetime signature 
(+—)• 


2 A solution representing a distorted Kerr-Newman black hole in 
an external static gravitational field was constructed by Breton, 
Garcia, Manko, and Denisova fl5t . 
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During each step we slowly bring the sources towards 
the black hole in such a way that the space-time axial 
symmetry is preserved. 

In such a case, the generated gravitational and elec¬ 
tromagnetic waves are weak and do not carry angular 
momentum. These waves accompany a transition from 
one value of the distortion field to the next that is in¬ 
finitesimally nearby. Such a transition can be considered 
as a perturbation of the initial state defined by the met¬ 
ric g£l and the electromagnetic 4-vector potential A o' 1 . 

( f) 

The result of one transition step is the final state t/W 
( f) 

and Aa ; which can be expressed as follows: 

9ap — 9al + ^^a/3 + 0(A 2 , Xg, p 2 ), (1) 

A^ = A^+gB a +0( A 2 ,Ap,p 2 ), (2) 

where at the black hole horizon h a p = 0{g a p ) and 
g a PB a Bp = 0(1), and where A, p <C 1 are dimension¬ 
less parameters defining the strength of the gravitational 
and electromagnetic perturbations, respectively. 

The perturbation during one transition step grows 
from zero to the final value corresponding to the new 
distortion field during a finite time interval. After the 
perturbation reaches its final value, the gravitational and 
electromagnetic waves decay as an inverse power of the 

advanced Eddington coordinate v 0,113 , and the black 

hole settles down to another state defined by new values 
of the parameters Di. Note that according to the na¬ 
ture of the perturbation, the black hole’s electric charge 
Q and angular momentum J remain fixed. We are as¬ 
suming classical positive energy conditions that do not 
allow the black hole event horizon area A + to decrease, 
so for transitions that are adiabatic, A + must remain 
unchanged. 

We shall show now that a sufficiently-slow quasi¬ 
stationary transition is indeed adiabatic, i.e. A + is 
an adiabatic invariant with respect to the perturba¬ 
tion. In order to do so, we shall follow the approach 
given by Hawking and Hartle 01 ■ Namely, we consider 
the Hawking-Hartle null tetrad (l a , n a , m a , fh a ) which is 
well-behaved on the future EH (see, e.g., 0)- At the 
horizon, the null vector l a = dx a /dt is a null geodesic 
generator of the null EH hypersurface. Here T is a non- 
affine parameter along the null geodesic generators of the 
EH, which we shall choose as one of the spacetime coor¬ 
dinates, so that l a = <5“. 

At the fixed t coordinate the EH surface is topologi¬ 
cally a sphere on which the orbits of the Killing vector are 
closed circles (except for two fixed points on the symme¬ 
try axes). We choose (f> to be the azimuthal Killing coor¬ 
dinate. Then, one can always choose the other coordinate 
9 on the sphere to be orthogonal to the <f> coordinate. 

We choose rasa coordinate which is constant (x = 0) 
on the EH and orthogonal to the 9 and <j> coordinates. 
The complex null vector m a and its complex conjugate 
ro“ for x = 0 lie on the EH. There is gauge freedom in 
their spatial and null rotations. We fix the null rotation 


gauge by imposing m 4 =0. To fix the spatial rotation we 
take mP real and mP imaginary. Given these null vectors, 
the null vector n a is then uniquely defined through the 
null tetrad orthogonality conditions in which the only 
nonzero dot products are l a n a = —m a fh a = 1. In the 
given coordinates we have the null tetrad on the EH, 

l a = [1,0,0,0], n a = [n*, n x , 0,0], 

m a = [0,0, e a , ie b ], fh a = [0,0, e a , — ie b ], (3) 

where the vector components are real functions of the 
coordinates t and 9 on the horizon. 

In this null tetrad the first two of the Newman-Penrose 
equations (see 0, 0|) take the following form at the 
black hole EH: 


p 2 + aa + 2 ep + T 00 , 

(4) 

2per + 2e<r + T 0 , 

(5) 


where p = lp- tOL m a mP is real (because l a is hypersurface- 
orthogonal) and measures the convergence of the null 
geodesic generators, and a = l a - i pm a m& is the shear, 
which is complex. For the tetrad ([3]), e is real, and be¬ 
cause l a is a null geodesic generator, one has l a lp- a = 
2 dp, which fixes the scaling of l a on the EH in terms of 
e. For the null tetrad ([3]) these spin coefficients read 

1 1 n x t 

P= j( a ,t + b, t ) i cr = -(a t -b t ), e = ~^~ ( 6 ) 

The complex Weyl scalar To = C a p 1 5 l a m ,3 l 1 m s corre¬ 
sponds to an ingoing transverse gravitational wave. The 
transverse directions of the wave defined by the vectors 
m a and m“ are along the black hole horizon surface, and 
the wave propagates in the direction defined by the ingo¬ 
ing null vector n a . Using the expressions © and © we 
derive 

fo = - b,tt - a 2 + b 2 t ) - e(a t - 6, t ) ■ (7) 

An ingoing electromagnetic wave is defined by the 
complex Ricci tensor component $oo = 2<^o<(>o = 

—(1/2 )R a pl a lP, where = F a pl a m^ is the Maxwell 
scalar. In the radiation gauge A t = 0, 

4>o = Ag^e a + iA^^eP , Too = 2 (Ag t e 2a + A^e 26 ). (8) 

Given the quantities above, one can calculate the rate 
of change of the EH area (cf. 0), 

< 9 » 

where here and henceforth an integral over dA + denotes 
an integral over the closed two-dimensional surface of the 
EH at time t. Using the transition expressions m-m we 
can solve the Newman-Penrose equations perturbatively 
and find p, in order to calculate the EH area change. 
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To construct the gravitational and electromagnetic 
perturbative expansions we observe that the gravita¬ 
tional perturbation $0 is of the first order in A and the 
electromagnetic perturbation is of the first order in p. 
The gravitational perturbation induces through its non¬ 
linear interaction with the background electromagnetic 
field the term in cf >o which is of the first order in A. As a 
result, the corresponding Ricci tensor component $oo is 
a quadratic expression in A and p. These terms induce 
terms of the corresponding order in d> 0 , and so we get 
the following power series expansions: 

d'o = A'kg 1 ’ 05 + 0(A 2 , A/x, /x 2 ), (10) 

4>o = ak/>o° ^ + A^o ’ + 0(p 2 ,Xp, A 2 ), (11) 

$00 = p 2 ^ + ^oo 1] + A 2 *& 0) 

+ 0(/x 3 ,A/x 2 ,A 2 /x,A 3 ). (12) 

The corresponding expansions of the spin coefficients 
can be deduced from the Newman-Penrose equations. 
According to the order of $o and $oo, the leading terms 
of the convergence are quadratic in A and p, while the 
leading term of the shear is linear in A, 

p = A 2 p( 2,0 i + A/xp^’i) + /x 2 p(0’ 2) 

+ 0(p 3 , \p 2 ,x 2 p,x 3 ), (13) 

cr = Aer (1,( h -|_ 0(A 2 , Xp, p 2 ) , (14) 

e = e (0) + Ae (1 ’ 0) + /xe ((U) + 0(A 2 , A/z, p 2 ) . (15) 

One can show that according to the equation (l a l a ) ; p = 
—2nlp defining the surface gravity n when the metric is 
stationary, we have 2e(°) = n. 

Substituting the expansions m-m into the 
Newman-Penrose equations and solving them we derive 
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a 

1 

r 

P+- 

o^- 

0 

(16) 


pOO 

J e-^'-V (\<T {1 ’ 0) (t')\ 2 + <f> { oo°Xt' 

)) dt', 



(17) 

■ (1 ’ 1} w = - 

pOO 

J e-^'-V^XtXdt 1 , 

(18) 

t°' 2 Xt) = - 

poo 

J e-^’-^^XtXdt'. 

(19) 


Here the convergence Pg 2 ’°^ is due to the primary gravi¬ 
tational and induced electromagnetic waves, the conver¬ 
gence Pq 1 ' 1 ^ is due to the primary and induced electro¬ 
magnetic waves, and the convergence p(°’ 2 ) is due to the 
primary electromagnetic wave. These convergences give 
the corresponding rate of change in the EH surface area. 

According to our model, for one transition step the 
gravitational and electromagnetic perturbations begin at 
t = ti and end at t = tf. Note that for t < ti, a and p 
are non-zero but exponentially small. This is because the 
solution HUD-HUD i s defined by the entire future history 


of the black hole, what is a manifestation of the teleolog¬ 
ical nature of the EH. As a result, for the effect of one 
transition step we derive 

^ = —2 Jp^\t)dA + , (20) 

/ +oo n 

dt p( net \t) dA+ , (21) 

-OO J 

p {net \t) = AV 2 > 0) (t) + A pp (hl \t) + mV°’ 2) (1) -(22) 

The total quasi-stationary transition consists of N such 
steps. 

Let A g = \gffl\gap - g^pl and a em = 
^J\g^(A a — Aa^)(Ap — A^)| be the total dimension¬ 
less perturbations of the original metric and of the 

electromagnetic potential A& representing an undis¬ 
torted Kerr-Newman black hole which in a quasi- 
stationary transition is brought to the distorted metric 
g a p and the corresponding electromagnetic potential A a . 
Then one gets A ~ A G /iV, p ~ Aem/N , and the total 
time of the quasi-stationary transition is of the order of 
NSt, where 8t = tf — ti- 

To get an estimate of the total area change A A + we 
model the perturbations H/g and $ 0 o by rectangular im¬ 
pulses. One can define an upper bound crude estimate 
for the amplitudes of the perturbations $0 and $00 ■ It 
follows from Eqs. HZD and ([ED, and from the order of the 
surface gravity k = 2e ^, that on the EH the order of the 
amplitude of $0 is 1 /(St) 2 for a fast transition, nSt -C 1, 
and is k/S t for a slow transition, nSt 1, while the 
order of magnitude of $00 is 1 /{St) 2 , regardless of the 
transition rate. 

One can check that in order to have the first terms 
on the right hand sides of the Newman-Penrose equa¬ 
tions negligible in comparison with the following ones, 
as needed in the expansions USD-HUD, one should have 
A -C nSt and p <C ( nSt ) 3 / 2 for nSt 1 and X, p -C nSt 
for nSt 1 (in which case we have already assumed 
the stronger condition A,/i<Cl). Substituting this per¬ 
turbation into the expressions (fT6l) - ([HJl) and HUD, we can 
estimate the order of the total relative change in the hori¬ 
zon area during the total quasi-stationary transition for 
small and large values of St compared to 1 /k, 

~ (A g/( k dt) + A g A em + A| m ) > 

(23) 

U<5i»l ~ (A g + A G A em + A%m) ■ (24) 

Thus, for A + to be an adiabatic invariant, i.e., for 
A A + /A + <C A g with respect to the quasi-stationary 
transition, we must have 

> (A . G /(/c<$t) + A em + A ^ m/Ag) )^’ (25) 
A^| K( 5t»i ^ ( A g + A em + A \ m ! A g )-77 ■ (2d) 
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For given values of Ac and A em, these expressions show 
the order of the minimum number of steps for the quasi¬ 
stationary transition for A + to be an adiabatic invari¬ 
ant. Note that because the electromagnetic perturbation 
would be expected to produce a gravitational distortion 
at least Ac ~ A EM , we expect that A em /Aq < 1 . 

The next step is to consider the fate of the black hole’s 
CH. As it was demonstrated by Poisson and Israel [22 1 . 
as well as by Ori [TEH , the infinite blueshift of a radia¬ 
tive tail produced at the CH results in the formation of 
a weak null curvature singularity, the mass inflation null 
singularity , at the ingoing part of the CH. In a follow¬ 
ing investigation by Marolf and Ori E3 , it was shown 
that for late-infall-time observers (the observers whose 
Eddington’s advanced null coordinate v e h at which they 
cross the EH, with v e h = 0 at the formation of the EH, 
obeys nv e h 1), the space-time in the black hole interior 
may be described by an essentially unperturbed station¬ 
ary space-time up to very near the CH. At the outgoing 
portion of the CH the linear gravitational metric per¬ 
turbations decouple into a superposition of outgoing 
and ingoing components and for w > 1 the outgoing 
component decays as an inverse power of v. 

The late-infall-time observers who arrive at the out¬ 
going portion of the CH will encounter a gravitational 
shock-wave singularity: a finite jump in the metric 
Ah,^ ~ 1, which is the shock wave amplitude, within 
an effectively vanishing observer’s proper time (the shock 
wave width) oc e~ KVeh <C 1, where k is the CH surface 
gravity. All the time-dependent structure of the space- 
time gets compressed into the shock wave. The shock 
wave singularity is more violent than the mass infla¬ 
tion null singularity , and it is detected by both freely 
falling and accelerated late-infall-time observers. The 
late-infall-time observers who arrive at the ingoing por¬ 
tion of the CH will encounter the mass inflation null sin¬ 
gularity which decays as an inverse power of v e h and van¬ 
ishes in the late-infall limit, v e h —> + 00 . 

This analysis was done for the Reissner-Nordstrom 
space-time. However, in the papers above arguments 
were given that the same situation will take place in the 
case of a rotating Kerr-Newman black hole (see discus¬ 
sions in El,S3)- This result suggests to us to consider, 
instead of the singular CH, a regular stretched Cauchy 
horizon , which is located very near the CH, which keeps 
track of the perturbation. A replacement of the CH with 
the stretched horizon is done by identifying events on the 
stretched horizon, e.g. an ingoing beam of light, with the 
corresponding events detected by freely falling observers 
very shortly before they enter the CH. The stretched hori¬ 
zon area A_ is approximately equal to the area of the 
would-be regular CH. The difference between these areas 
decreases exponentially with the advanced time v e h after 
the perturbation and vanishes in the late-infall limit. 

Here we shall assume that the results of Marolf and 
Ori E2 are valid for a distorted black hole as well as 
for the undistorted black holes they studied. Then one 
can define and solve the Newman-Penrose equations ([3])- 


([5]) at the stretched Cauchy horizon and take the limit 
to approach the real Cauchy horizon from the regular 
interior region and use the equation for the area rate |3|). 
As a result, A_ is an adiabatic invariant as well. Thus, at 
the end we have both the area A_ of the CH and the area 
A + of the EH remaining adiabatically constant, so that 
the relation A+A_ = (8nJ) 2 + (47 tQ 2 ) 2 , which is true for 
an undistorted Kerr-Newman black hole, remains valid 
for distorted black holes as well. One can also conclude 
that not only the area product but any function of both 
A + and A_ remains the same. 

The horizon area product for higher-dimensional black 
objects was studied as well. For example, Castro and 
Rodriguez [24| showed that product of horizon areas is 
independent of the mass and topology of the horizons for 
all known 5-dimensional asymptotically flat black rings 
and black strings. However, as it was illustrated by Visser 
[25l | . for different multi-horizon black hole solutions not 
all area products are mass independent. In the work 
by Cvetic, Gibbons, and Pope jl|, explicit results for the 
product of all horizon areas for general multi-charge black 
holes were given in terms of quantized charges, angular 
momenta, and the cosmological constant @ 

All these solutions are exact, without influence of ex¬ 
ternal matter and fields on the black objects. Thus, it 
is important to know whether these area products are 
generic as they are in the case of the Kerr-Newman black 
hole. Therefore, it is natural to ask whether our heuris¬ 
tic argument can be generalized to higher-dimensional 
black objects. We can follow the same scenario and show 
that, assuming a black object is stable against an external 
perturbation due to a (higher-dimensional) weak gravi¬ 
tational wave generated by motion of distorting matter 
which preserves some symmetry (an analogue of the axial 
symmetry) such that the black object’s angular momenta 
(as well as charges) are preserved, the outer horizon area 
does not change due to an adiabatic perturbation. As 
far as the inner horizons, some arguments were given in 
E2 that the late-time limit of black hole interiors can 
be carried over to generic stable black holes in any di¬ 
mension. This implies that we can repeat our scenario 
by introducing the stretched inner horizon(s), as we did 
in the case of the Kerr-Newman black hole, and conclude 
that the inner horizon area(s) do not change either. 

Thus, we can formulate the following 
Conjecture: Assume one has a d-dimensional (d > 4), 
stationary, multi-horizon black-object solution with the 
horizon areas {Ai,i = 1, ...,n} ; a set of charges Qj, and 
angular momenta {Jk,k = 1, ...,[(d — 1) /2]}, and possibly 
with a cosmological constant A, such that the following 


3 Note that it was argued by Faraoni and Zambrano Moreno [26ft 
that quantization rules for the horizon areas of stationary black 
holes are misleading, for they do not correspond to realistic time- 
dependent situations. 
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relation holds: 

n 

n^ = m,4,A- i/2 ), 

i=1 

where f is some polynomial function. Then, assuming 
that the solution is stable against external gravitational 
perturbations, this relation will hold for any adiabatically 
distorted solution with the same values of the charges and 
angular momenta. 

As an example illustrating this conjecture, we can con¬ 
sider a five-dimensional, static, electrically charged black 
hole distorted by an external, static and electrically neu¬ 
tral distribution of matter. Then the product of the ar¬ 
eas of the outer and inner horizons is proportional to the 
cube of the electric charge j2?j, Hl| ■ 

One of the possible issues related to this conjecture is 
a consideration of the inner horizons of a multi-horizon 
solution which are located behind the outermost inner 
one (assuming that the solution has nested inner horizon 
topology.) Due to the gravitational perturbation of the 
outermost inner horizon, the mass inflation phenomenon 
and formation of a gravitational shock wave singularity 
will seal off the interior region. As a result, only two areas 
(of the outer horizon and the outermost inner horizon) 
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